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Abstract 

Particle multiplicities and ratios in the microcanonical ensemble of relativistic gases near production 
thresholds are studied. It is shown that the ratio of heavy to light particle multiplicity may be enhanced 
in comparison to its thermodynamic limit. 
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I. INTRODUCTION 


The statistical model of hadron production in high energy collisions of both elementary and 
heavy ion collisions proved to be remarkably succesful in reproducing the multiplicities of dif¬ 
ferent hadron species. This hnding has raised the question of the meaning of thermodynamical 


quantities in these systems and triggered a still ongoing debate [1|, 

In the case of relativistic heavy ion collisions, due the large multiplicities involved, the use 
of the grand canonical ensemble (GCE) is sufficient and the model predicts that hadron yields 
should be linearly dependent on the volume at hadronization. On the other hand, in elementary 


collisions, the canonical ensemble (CE) [l| 


enforcing the exact conservation of charges 


must be used owing to the relatively low multiplicities. At lower multiplicities, even the micro- 
canonical ensemble (MCE), where energy-momentum conservation is enforced, is needed to 
correctly describe particle abundances. The transition from MCE to CE has been studied in 


>criDe pa 

0 , 0 - 


detail in refs. 13, ISl. The equivalence between the three ensembles, with respect to particle 
multiplicities, is recovered only in the thermodynamic limit, when U —> cx). 

The charged conservation laws in the CE imply a reduction of the mean multiplicity of heavy 
charged hadrons in comparison to the corresponding multiplicity calculated in the CCE. This 
effect is usually defined as canonical suppression and was particularly studied for the production 


of strange 0 and open charm hadrons as well as anti-baryons 

The canonical suppression is significant for the charged particles with a mean multiplicity 
smaller than 1. Naively one may expect that an introduction of the energy conservation law by 
use of the MCE should lead to the gradually increasing suppression of the multiplicity of heavy 
hadrons with decreasing energy of the system. Contrary to this expectation, in this work we 
demonstrate that the energy and momentum conservation laws lead to rich structures in the 
dependence of hadron multiplicities on the system energy. 

The paper is organized as follows: in Section II we study this effect by means of a simple 
system of massless and heavy particles, where an analytical description is possible. In Sect. 
Ill, we study the full ideal hadron-resonance gas by means of numerical methods described in 


a 


refs. 


Conclusions are summarized in Sect. IV. 
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II. THE ANALYTICAL MODEL 


In this section we consider the system of noninteracting Boltzmann particles with zero mass, 
both neutral (with degeneracy factor g) and charged (with degeneracy factor g±). The system 
also includes neutral heavy particles with mass M and degeneracy factor G. The nonrelativistic 
approximation is used for heavy particles. 


A. Grand Canonical Ensemble 


In the GCE independent variables are volume V, temperature T, and chemical potential /i. 
The average number of massless neutral particles N, massless charged particles, iV+ and N_ , 
and heavy neutral particles, n, are: 


(N) 

(n) 


GCE — 


gVT^ 


vr^ 


, {N±)gce — 


9±VT^ 


GCE 


= GV 




TT^ 


exp(±^). 


exp 


M\ 

~Y 


( 1 ) 


The system energy reads: 


{E)gce = — ^T{N)gce + 3T{N^)gge + 3T{N_)gce + ( -^T + M ) {n)GCE- (2) 


Of course, the nonrelativistic approximation used for massive particles holds for M ^ T. 


B. Micro-Canonical Ensemble 


In the MCE independent variables are volume V, energy E, and net charge Q. The MCE 
partition function fl]s[{E, V) in the system of N neutral massless particles equals to ^ : 

N noo noo / ^ \ -i / _t;"\ N 


= T(|h) ^ pi,dp^s(E-^^p^ ^ 


(3A- 1)! ■ 

( 3 ) 

Let us consider the case of charged massless particles with total zero net charge, Q = 0. This 
implies /xq = 0 in the CCE and thus {N+)gce = {N-)gce- The MCE partition function for 


iV_|_ = N_ charged massless particles reads 


CatJE, I/) = 


/LLCl 

i: 




2N± 


j^6N± 


-1 


NY V y (6A± - 1)! ’ 


( 4 ) 
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where N± = N+ = N_. 

The MCE partition function for n heavy nonrelativistic particles can be calculated analyti¬ 
cally: 




pldpi... pldpns 


E -^(m + 
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2M 
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Ti r 
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{E - nM) 2 
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(5) 


n\ \(27r)^/^/ r (^) 

where the Euler gamma function r(a;) has a simple form for integer k and halfinteger k + 1/2 
arguments: 

1 ■ 3 ■ . .. ■ (2fc - 1) 


m = (fc- 1 )!, r(fc + -) = 


TT . 


Using Eqs. (ElED one can calculate the partition function for N massless and n heavy particles: 


U) = 


noo poo 

/ dE, / dE^ n^{E,, V) V)5[E - E, - E^] 

'o Jo 

N 

/ f -I / \ i\ /I o o 
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1 1 fgV 
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, ^ 3 n 
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( 6 ) 


Note that the MCE partition functions Qn^E, V) 0 and fln{E, V) ® are defined for nonzero 
particle numbers, > 1 and n > 1, because of an exact energy conservation. However, the 
MCE partition functions Q]sf^n{E, V) (jHI) requires only -|-n > 1, so that either N or n can be 
equal to zero. 

Finally, using Eqs. dllSD one hnds the partition function for N massless neutral, = N_ 
massless charged, and n heavy particles: 


nN,N^AE, V) = dE, dE 2 QnAEi, V) nNAE2, v) s[e -e,- e^] 
Jo Jo 
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TV! (iV±; 
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X 


r(3A^ 6A^± |n) 
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(7) 


The MCE partition function r2Ar,Af±,n(-E, V) (|7j) is evidently transformed into fl^AE, V) (P) 
for g±. = 0. For g = 0 the V) from Eq. (|7j) is transformed into the partition function 

for = N_ massless charged and n heavy particles: 


^m^AE, V) 


1 1 fg±vA^^ / GV 

(Ayp^dV^y V(27r)3/2 


■71 3n 

M 2 


T{QN± + In) 


{E - nM) 
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( 8 ) 
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The MCE partition function can also be calculated for the case of non-zero system net charge 
Q in two steps: 


V) 




(El, E)fiiv_(E 2 , V)6[N+-N_ 


Q]5[E-E^-E2], 

(9) 


nN+,N.,n{E,V) = / dE, dE 2 f]iv+,iv_(Ei, E)<5 [E-Ei-E2]. (10) 

Jo Jo 

Finally, the partition function for positively charged massless particles and for n neutral 
heavy particles equals to: 


C, Q) = 


1 1 1 (g±v\ 

N+\ {N+-Q)\ n! V vr2 ) 
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r( 6 Ar+ - 3g + |n) 
and for negatively charged massless particles: 

^N.,n{E, V, Q) = 


(E - , 


1 (9±V\ 
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( 11 ) 


( 12 ) 


These formulas transform into VLM±,n{E, V) ((HI) for Q = 0. 

Now we calculate the mean particle multiplicities in the MCE and compare them to the 
GCE multiplicities m- This comparison is performed at the same energy density which is set 
to be the energy density calculated in the GCE for T = 160 MeV. In the MCE the volume 
V equals to that in the GCE, and the MCE energy E equals to GCE average energy Egce 
0 . Let us start with the system of massless charged and heavy neutral particles with zero net 
charge. In this case all MCE averages are calculated as: 


(• • ■)mce 


1 

G(E, E; Q 


0, M) 


...G^^,„(E, E), 

N±,n=0 


(13) 


where Gv±,n (E, V) is given by Eq. 0, Vl{E,V]Q = 0, M) = Yl'N±,n=o^N±,n{E, V), and 
r2o,o(-£^; ^) = 0, Umax = [^] due to the exact energy conservation. The degeneracy factors 
are fixed to be 5 ^^ = G = 1. In Fig. Qthe ratios of average particle numbers in the MCE to 
those in the GCE are presented as functions of the system energy E = {E)gce- We remind 
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that V = E/e{T). Circles and triangles represent Monte-Carlo (MC) calcnlations for the same 
parameters (see Section III for details of the nsed procednre). These calcnlations are in an 
agreement with the analytical resnlts. Fnrther we also show the MC resnlts obtained with the 
additional reqnirement of the exact momentnm conservation (P = 0) for the whole system. For 
the heavy particles with mass M = 3.1 GeV one observes (see Fig. the MCE snppression of 
the heavy particle mnltiplicity. This snppression is strong at energies E comparable to heavy 
particle mass M, and it is still signihcant aX E = 100 GeV S> M. 

In Fig. El (left) the ratios {N±)mce/{N±)gce and (n) mce/( n)GCE are shown for a 
smaller valne of the heavy particle mass, M = 0.7 GeV. One observes a fast drop of the 
{N±)mceI{N±)gce and strong “oscillations” of the (n) mce/{ n)GCE ratio. The MGE enhance¬ 
ment effect {{n)MCE/{n)GCE > 1) increases strongly with decreasing mass M of the heavy 
particle. Note that the non-relativistic treatment of the heavy particles is still approximately 
valid at M = 0.7 GeV and T = 160 MeV. 

Let us consider the massless charged and heavy neutral particles with the net charge Q = 2. 
In this case all MGE averages are also calculated by means of Eq. (IT!?|l . Besides, {N+)mce and 
(7V_) MCE are calculated independently using Eqs. mm . The corresponding GGE values now 
contain exp(±/i/T), and the equation {N+)gce — {N-)gce = Q should be solved: 


9±VT^ 




Q. 


(14) 


where the volume is given by Eq. (El) as {E)gce/£{T). In order to make a simple estimate, we 
neglect a small contribution of heavy particles to the energy density and consequently Eq. dn 
transforms into: 


— tanh 
3T 



(15) 


This equation has solutions for E > 3T Q. The results for the MGE to GGE particle ratios are 
presented in Fig. El (right). Negatively charged particles are suppressed because of the exact 
charge conservation, and the hrst maximum seen in Fig. El (left) for heavy neutral particles 
disappears. 
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C. Canonical Ensemble 


Let us consider the CE system with zero total net charge, Q = 0. A standard way is 
to coi^are the CE and GCE results at the same volume V and temperature T (see e.g. 
1^1 IM 1^)- One observes here the reduction of charged particle multiplicities with decreasing 
system volume in comparison to ones calculated in the thermodynamical limit. The effect is 
called canonical suppression. Note that the energy density in the CE, Ece, is also suppressed 
in comparison with that in the GCE at the same temperature. This is due to the suppression 
of charged particle multiplicities. Our aim is to compare the particle multiplicities in different 
ensembles at the same energy density. Thus, we perform the comparison at the same volume 
and energy in all ensembles, so that the energy densities are also equal to each other, Ece = 
£mce = £gce = E/V = const. Of course, under this condition the mean multiplicities and 
temperatures calculated in the CE and the GCE are equal in the thermodynamic limit, V oo. 
The results are, however, different for small systems. When the volume of the system in the 
CE decreases, but the energy density is kept hxed, the temperature increases which leads to 
an increase of the heavy neutral particle multiplicity. The CE temperature, T*, is obtained as 
a solution of the equation: 


^ QT* + = y 

where {z = gVT*^ 


^GCE — const , 


(16) 


{N±)ce — z 


h{2z) 

Io{2z) 


{E)ge — GV 


MT* 

271 


3/2 


exp — 


M 


(17) 


A comparison of Eq. (HZD and Eq. m demonstrates the CE suppression of charged multiplicities. 
This leads to T* > T due to Eq. (ITT)|l hence to {n)cE > {n)GGE , as the equations for neutral 
particle multiplicities have the same form in the GCE and the CE. The results are shown 
in Fig. El Massless particles are suppressed in the CE, Fig. El (left), since the increase of 
temperature is insufficient to overcome the CE suppression. A smooth increase of the ratio 
{'ci)ceI{'cl)qce is seen in Fig. El (right). Thus, a fast increase just above threshold and the 
following “oscillations” of the {n) mge/{E)gce ratio are due to energy-momentum conservation 
laws. 
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D. Grand Micro-Canonical Ensemble 


We turn now to the discussion of the system with neutral particles only, both massless and 
heavy ones. When the energy is exactly conserved we call this ensemble the grand micro- 
canonical ensemble (GMCE) in order to distinguish it from the MCE where both the energy 
and the charge are strictly conserved. The calculation within the GMCE allows us to check 
whether the MCE enhancement is only related to the charge conservation or it is due to the 
energy conservation. We set g± = ^ and g = G = 1, and in this case all GMCE averages are 
calculated as: 

-y OO 

“ fi(£, V;M) ^ 

^ ’ ’ N,n=0 

where ^1{E,V;M) = YIn n=o^ri,NiE, V) is the total MCE partition function, where 
f2o,o(-S, E) = 0 and Umax = due to the exact energy conservation. The result is shown in 
Fig.m where one can see that the heavy particle enhancement near the threshold is also present 
in the neutral system. The amplitude of a drop in {N)gmce/{N)gce and the “oscillations” of 
{n) GMCE!{n )gce becomes smaller than in the MCE with charged massless particles, but they 
are still present. The Monte-Carlo (MC) calculations with the total momentum conservation 
show a similar behavior: a decrease of the ratio for the massless particles and an increase for 
the heavy particles. The momentum conservation leads to an increase of both ratios, it also 
smears and shifts the peak for heavy particles. In Fig. El we present the energy dependence of 
the particle multiplicities calculated in the GMCE and the GCE. 


III. THE HADRON-RESONANCE GAS MODEL 


This section starts with a brief presentation of the numerical methods used to calculate mean 


hadron mult ip 


be found in 


icities in the MCE for the ideal hadron-resonance gas. The full description can 
Further on, we show and discuss the results on the energy dependence of 


hadron yields near the threshold. 





A. The Procedure 


The MCE partition function is defined as the sum over all multi-hadronic states \ hy) localized 
within the volume V of the system and constrained with the four-momentum and the abelian 
(i.e. additive) charge conservation: 


a = - P^)iQ,Qjhv }. (19) 

hy 


where Q = (Qi, • • • ,Qm) is a vector of M integer abelian charges (electric, baryon number, 
strangeness etc.), P is the overall four-momentum of the system and Pop, Qop the relevant 
operators. Provided that relativistic quantum held effects are neglected and the volume of the 
system is large enough to allow the approximation of hnite-volume Fourier integrals with Dirac 
deltas, it can be proved that the micro-canonical partition function D can be written as a 
multiple integral: 


'■+7r 


Q = 


(27r)4+^ 


y / d^0e 


iQ-0 


exp 


[ 5 : 


(2J, + i)r 

(2,r)l j 


d=^p log(l ± (^0) 


where (y is the vector of the abelian charges for the hadron species, Jj its spin and pj = 
+ P^, p), where rrij is the mass of the hadron species; the upper sign applies to 
fermions, the lower to bosons. The integral (PU|) is more easily calculable in the rest frame of the 
system where P = 0. Unfortunately, an analytical solution with no charge constraint is known 
only in two limiting cases treated in previous sections: non-relativistic and ultra-relativistic 
(i.e. with all particle masses set to zero). The full relativistic case has been attacked with 


several kinds of expansions but none of them proved to be fully satisfactory as the achieved 
accuracy in the estimation of different kinds of averages could vary from some percent to a 
factor 10. Therefore, a numerical integration of Eq. (PU|) is needed. The most suitable method 
is to decompose D into the sum of the phase space volumes with hxed particle multiplicities 
for each species: 

^ jv.qi > (21) 

{Nj} 


{iVj } being a vector of K integer numbers (A^i,..., Nk), i.e. the multiplicities of all of the K 
hadronic species. 

Here we use the approximated expression in Boltzmann statistic of the phase space volume 
for the channel {Nj} in order to make a qualitative study on micro-canonical effects in a 
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coherent way with the analysis in the analytical model made in the previous part of the paper. 
For the same reason, we will keep hxed masses for hadron resonances neglecting the Breit 
Wigner broadening effect. Both of these simplihcations, do not affect the results signihcantly. 


The phase space volume for hxed multiplicities in Boltzmann statistic reads 


Q: 


^{Nj} = 


n 




f d^pi...d^p7v S^{P-^pi), 

i=i 


( 22 ) 


where N = Y,j Nj. 

We are mainly interested in the calculation of quantities relevant to particle multiplicities, not 
to their momenta, namely their kinematical state. The average of an observable O depending 
on particle multiplicities in the micro-canonical ensemble can then be written as: 

NjCij 


{ 0 ) = 


(23) 


where 0({A^j}) = W for the average multiplicity of the fc—th hadron species (W)- Altogether, 
what we need to calculate in order to evaluate an average (ESI) are integrals like (El- 

In order to be able to effectively calculate fl{Nj} for any channel, a brute force option 
is to do it for all of them. However, this method is not appropriate for a system like the 
hadron gas, as the actual number of channels is very large. Indeed, with 265 light-havored 
hadrons and resonances (those included in the latest Particle Data Book issue llj), the number 
of channels allowed by the energy-momentum conservation is enormous and increases almost 
exponentially with the cluster mass, which requires an unacceptably large computing time. 
Therefore, the calculation of the phase space volume of all allowed channels is only possible 
for very light clusters, in practice lighter than ~ 2 GeV. Hence, if a method based on the 
exhaustive exploration of the channel space is not affordable, one has to resort to Monte-Carlo 
methods, whereby the channel space is randomly sampled. 

An estimate of the average ESI) can be made by means of the so-called importance sampling 
method. The idea of this method is to sample the channel space (i.e. the set of integers Nj, 
one for each hadron species) not uniformly, but according to an auxiliary distribution HiAr^.} 
which must be suitable to being sampled very efficiently to keep computing time low, and, 
at the same time, as similar as possible to the distribution The latter requirement is 

dictated by the fact that ^{Nj} is sizable over a very small portion of the whole channel space. 
Thus, if random conhgurations were generated uniformly, for almost all of them D{Ar^.} would 
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have a negligible value, and a huge number of samples would be required to achieve a good 
accuracy. On the other hand, if samples are drawn according to a distribution similar to 
little time is wasted to explore unimportant regions and the estimation of the average ()2d|l is 
more accurate. A crucial requirement for is not to be vanishing or far smaller than ^{Nj} 

anywhere in its domain in order not to exclude some good regions from being sampled, thereby 
biasing the calculated averages in a hnite statistics calculation. 

A Monte-Carlo estimate of (O) is: 


{ 0 ) = 


l^k=l 




n 

n 


(fc) 

wi 

(fc) 

{Nj} 


lA'} 

2^k=i (k) 


(24) 


where are samples of the channel space extracted according to the distribution If and 

fulhlling the charge constraint Q = Njqj. 

Provided that Ns is large enough so that the distributions of both numerator and denomi¬ 
nator in Eq. are Gaussians (hence the conditions of validity of the central limit theorem 
are met), the statistical error cr(o) on the average (O) can be estimated as: 


2 

^(O) 




-n 


O 


Q2 

TT^ 


+ (0)'E; 
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'q2 

^Va, 


2(0)En O 


Q2 

^Va 

^Va, 


(26) 


where En stands for the expectation value relevant to the fl distribution. 

Here we dehne Hijv^.} as the product of K (as many as particle species) Poisson distributions: 


K Nj 

= (26) 
j=i y 

which will be henceforth referred to as the multi-Poisson distribution or MPD, enforcing as 
mean values the mean hadronic multiplicities Uj calculated in the GCE with volume and mean 
energy equal to the volume and mass of the system: 


A = 


(2J,- + 1)1/ 

27r2 


m^TKs 


?)n+. 


(27) 


where V is the volume of the system, T is the temperature and A* the fugacity corresponding to 
the charge Qi. Temperature and fugacities are determined by enforcing the GGE mean energy 
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and charges to be equal to the actual energy E and charges Q of the system: 



3 * 




(28) 


3 * 


with 



(29) 


The distribution (EHl) can indeed be sampled very efficiently and is the actual multi-species 
multiplicity distribution in the GCE in the limit of Boltzmann statistics. 

B. Threshold Effects in the Hadron-Resonance Gas 

We have calculated the ratios of mean multiplicities of several hadron species in the MCE 
and the CE with respect to the GCE for an ideal hadron-resonance gas including all species 
up to a mass of about 1.8 GeV. Quantum statistic effects and resonance mass broadening have 
been turned off. In order to make a proper comparison between the different ensembles, energy 
(in the rest frame) and volume in both the CE and the GCE have been set to the same value 
as in the MCE. This implies that the energy density e has the same value in all ensembles as 
well. It has been hxed to 0.3895 GeV/fm^, corresponding to a temperature of 160 MeV in the 
GCE. As it has been discussed in Section II, in the CE, owing to the exact charge conservation, 
the energy density is not an intensive quantity but it also depends on the volume. This implies 
that for a hxed energy density, temperature varies as a function of total energy. In particular, 
for a completely neutral hadron gas which will be discussed henceforth temperature decreases 
as energy increases. This happens because at lower energies or volumes charge conservation 
suppresses more and more charged particle multiplicities and the system needs an increase in 
temperature to keep the energy density constant. 

The results of our calculations are shown in Figs. IHlfTT^ as a function of the total energy 
near the production threshold. In general, it can be seen from these plots that for an actual 
ideal hadron-resonance gas the (micro-)canonical enhancement for some particle species near 
the threshold (which was implicitly numerically observed in shows the same qualitative 
features as those of the analytical model. Nevertheless, the correlation effects between different 
particle species due to charge and the energy-momentum conservation play an important role 
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from the quantitative point of view and results in abrupt changes in the micro-canonical mean 
multiplicities as the total energy exceeds the threshold for the production of a new channel. 

In Fig. El we present the ratio between the average multiplicities of vr^-mesons in the MCE 
(CE) and in the GCE as closed circles (dashed line). Since pions are the lightest hadrons, we 
expect them to play the same role as massless particles in the analytical model. The ratios 
CE/GCE and MGE/GGE approach each other at energy of about 2 GeV and, thereafter, slowly 
converge to the GGE limit. For lower energies (less than ~ 2 GeV) the total particle multiplicity 
in the MGE is ~ 2, as the dominant channels are -|- 7r° and vr"'' -|- 7r“. Therefore, while in the 
GGE the mean multiplicity of 7r°-mesons is proportional to the energy, in the MGE it is almost 
constant and ~ 2/3. This explains why the ratio (7r°)MCE/(7r°)GCE has a nearly hyperbolic 
shape at small energies. 

On the other hand, the increase of the ratio (7r°)cE/(7r°)GCE as energy decreases is due to the 
higher temperature in the canonical ensemble for small systems with the same energy density 
(see Section irm and the discussion above). As the heavy (i.e. with m S> T) neutral particle 
multiplicity is proportional to exp[—m/T] according to Eq. (fT7|l . the ratio {nh)cE/{nh)GCE 
turns out to be proportional to: 

-rrih ( —r) j (30) 

\T{V,e) T{oo,e)J\ 

which grows exponentially with the particle mass. Indeed, in Figs. 17111111 the ratios for heavy 
neutral particles rj, p, /2 and J/'ip show the expected behavior, with a “canonical” enhance¬ 
ment larger for heavier particles. For J/t/, it is so large that it exceeds the “micro-canonical” 
enhancement, which, conversely, tends to diminish at larger masses. Indeed, for a sufficiently 
large mass, this enhancement disappears as it should have been expected from the result of the 
analytical model discussed in Section ITTl 

In fact, as in the case of the analytical model, we observe a strong enhancement of the 
average multiplicities in the MGE with respect to the GGE near the production threshold. 
This effect is larger for lighter particles: about a factor ~ 10 in the ratio MGE/GGE for the 
second lightest neutral meson rj and ~ 3.5 for J/-^ meson (see Figs. [7| and cnD- These are 
remarkable hgures, as the threshold peaks actually correspond to the “second” peaks observed 
in the analytical model (see Figs. El (left) and (right) as well as Fig. El (right)). The hrst peak 
(the largest one) is forbidden due to momentum conservation in the MGE that allows no less 
than two particles in the hnal state. A similar effect can be seen in the analytical model (Fig.Ej) 




CE 


(nh) 


GCE 


exp 
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for non-vanishing charge, because Q = 2 in the MCE also demands at least two particles in the 
final state. One can also see that the peaks in the full hadron gas have a gradual increase which 
is typical for the “second” peak (compare Figs. |21 (left) and (right)) in the analytical model. 

In full hadron gas 7r+ and 7r“-mesons play the role of the massless charged particles in 
the analytical model considered in Section II. The “second peak” of neutral heavy particle 
multiplicities is present similar to that in Figs. ISlandEl This “second peak” is stronger than 
in an analytical model, because the GCE energy density in the real hadron gas is a far larger 
for the same T = 160 MeV. This is due to much larger number of particle species. An exact 
momentum conservation makes the whole effect for the “second peak” even stronger. 

The exact energy-momentum conservation is responsible for a non-smooth variation of the 
mean multiplicities as a function of the energy of the system. In fact, the opening of a new 
channel at a total energy E implies a rapid change of the behavior of this function. In order to 
show this, it is appropriate to have a closer look at the MCE/GCE ratios for vr^, rj, p-mesons. 
In Fig. nria zoom thereof is shown on the energy range 0.6-1.5 GeV, where threshold energies 
for relevant channels are marked. It is evident that there is a full correspondence between 
them and the slope changes. The general pattern can be understood from the analysis of the 
analytical massless-heavy particle model. The MGE/GGE ratio of massless particles gradually 
decreases to the threshold of heavy particle production. At the threshold, massless particle 
yield drops and heavy ones are enhanced. Likewise, in the full hadron gas 7r°-mesons play the 
role of the lightest particles and always show a drop whenever a new channel is opened. All 
other particles behave like “heavy” ones near their threshold, but when a new heavier particle 
threshold opens up, their MGE/GGE ratio drops, just as the lightest particles would do. For 
instance, in the middle panel of Fig. the rj production sets in at the channel threshold 
and the curve raises until the energy m^rO -|- rrip is reached. Here the channel is opened and, 
just a little later, too. The production onset of these two new heavier particles downgrades 
7] to the role of a light particle and a decrease in the slope of the ratio {i])mce/{v)gce is implied. 
For little larger energy, the threshold for the r/r/ channel entails a new step up in the slope of 
the r] ratio. Then, for E > niri + rrip, the channels r]p and rjuj are open and, consequently, 
both T] and p ratios increase. Then, for energies larger than ~ 1.5 GeV, the number of allowed 
channels increases quickly and the ratio MGE/GGE shows a rocky pattern, which smoothes 
out thereafter. 
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Finally, we turn the quantities which might be, in general, related to the experimental results 
on hadron production in collisions of relativistic particles. Under simplifying assumptions that 
the hadron production is statistical and the system energy is proportional to the collision energy, 
hadronic hnal states can be identihed with different conhgurations in the MCE. Consequently 
the dependence of the ratios of the mean hadron multiplicities in the MCE on the system energy 
may be expected to be related to the collision energy dependence of these ratios. 

The dependence of (vr^) on the system energy calculated by use of the MCE and the GCE 
is presented in Fig. El There is a smooth and monotonic increase of (vr^) in the GCE. A very 
different behavior is seen in the MCE. For E < 3 GeV the pion yield reveals a non-monotonic, 
rocky dependence. In the very nearly of the threshold, (vr®) = 2. This is due to the small 
difference between the mass of 7r° and tt^ (~ 5 MeV) that prevent charged pions from being 
produced. For slightly higher energy, as the channel 7r+ -|- 7i~ is open too, (vr^) quickly decrease 
to ~ 2/3 as said before. 

In Fig. na the ratios {f])/{7i^), (p)/(7r°) and {J/'ip)/{7i^) are plotted as a function of the 
system energy. In the GCE the ratios are independent of the energy, whereas in the MCE a 
rocky pattern is predicted. Only above E > 5 GeV the ratios smoothly approach from above 
the GGE limit. 

In this work we made the assumption of a constant energy-density, that is, in every plot the 
volume is proportional to the energy of the system. Different hypothesis (for instance an equal 
constant temperature for GGE and GE and an energy-density for MGE: e{T) = E/{V)ce) 
would lead us to a completely different behavior of hadron multiplicities near the threshold 
with the same thermodynamical limit. A close comparison with the data could shed some light 
on which is the right critical quantity at hadronization. 


IV. SUMMARY 

We calculated and discussed the threshold behavior of particle multiplicities in relativistic 
gases. The micro-canonical formulations of two models were used. Firstly, analytical formulas 
were derived for the gas of massless and heavy particles. Second, the results for the hadron- 
resonance gas are obtained by use of the Monte-Garlo procedure. In both models the ratio 
of heavy to light particles near the threshold is enhanced in comparison to one in the grand 
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canonical limit. At low system energies {E < 3 GeV) a non-monotonic and rocky dependence of 
hadron multiplicities and their ratios on the system energy is obtained within hadron-resonance 
gas model. 

The problem whether these unexpected micro-canonical threshold effects can be seen in the 
experimental data on collisions of relativistic particles is left for a future study. 
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FIG. 1; The dependence of the ratios of the mean particle multiplicities calculated in the MCE to those 
the GCE on the total energy of the system. The gas with zero net charge Q = 0 of massless charged 
particles {g± = 1) and heavy neutral particles (M = 3.1 GeV and G = 1) is considered. The ratios for 
massless particles, {N±)mce/{N±)gce, (dashed line) and heavy neutral particles, {n) mce/ {n)GCE, 
(solid line) are plotted. The corresponding Monte-Carlo (MC) calculations are shown by circles and 
triangles. The calculations are performed for T = 160 MeV and E = {E)gce in the GCE and for 
E = s{T)V in the MCE. 
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FIG. 2: Left: The same as in Fig. Qbut for heavy particle mass M = 0.7 GeV. 

Right: The same as in Fig. ^but for heavy particle mass M = 0.7 GeV and the system net charge 

Q = 2. 
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FIG. 3; The dependence of the ratios of the mean particle multiplicities calculated in the MCE and 
the CE to those in the GCE on the total energy of the system. The gas with zero net charge Q = 0 
of massless charged particles {g± = 1) and heavy neutral particles (M = 0.7 GeV and G = 1) is 
considered. The ratios for massless particles, {N±)mce/{N±)gce (solid line) and {N±)ce/{I^±)gce 
(dashed line) are shown in the left plot whereas the ratios for heavy particles, {n±) mce/{ n±)GGE 
(solid line) and {n±)c e/{ n±) ggE (dashed line) are shown in the right plot. The calculations are 
performed for T = 160 MeV and E = {E)gge in the GCE and for E = s{T)V in the MCE and in the 
CE. 
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FIG. 4: The dependence of the ratios of the mean particle multiplicities calculated in the MCE to 
those in the GCE on the total energy of the system. The gas of massless neutral particles {g± = 1) 
and heavy neutral particles (M = 0.7, 1.3, 3.1 GeV and G = 1) is considered. Left: The ratios for 
massless particles, {N±)mce/{^±)gce calculated within the analytical model without momentum 
conservation are shown by solid (M = 0.7 GeV) and dashed (M = 3.1 GeV) lines. The corresponding 
Monte-Carlo results with momentum conservation are indicated by triangles and circles, respectively. 
Right: The ratios for heavy particles, {n±) mce/{ n±)GCE calculated within the analytical model 
without momentum conservation are shown by solid lines. The corresponding Monte-Carlo results 
with momentum conservation are indicated by triangles, squares and circles. The calculations are 
performed for T = 160 MeV and E = {E)gce in the GCE and for E = £{T)V in the MCE. Connecting 
lines between the MC-dots are drawn to guide the eyes. 


21 












FIG. 5: The dependence of the mean particle multiplicities calculated in the MCE, the GMCE and 
the GCE on the total energy of the system. The gas of massless neutral particles {g± = 1) and 
heavy neutral particles (M = 0.7 GeV and G = 1) is considered. Top panel shows the multiplicity of 
massless particles, {N±), calculated in the GCE (dashed line), in the GMCE (solid line) and in the 
MCE with momentum conservation (triangles). Bottom panel shows the corresponding results for the 
heavy particles. The calculations are performed for T = 160 MeV and E = {E)cce in the GCE and 
for E = e{T)V in the MCE. Connecting lines between the MC-dots are drawn to guide the eyes. 
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FIG. 6; The energy dependence of the ratio of the mean Tr^-meson multiplicity in the MCE and 
the CE to those in the GCE obtained within the hadron-resonance gas model. The calculations are 
performed for T = 160 MeV and E = {E)gce in the GCE and for E = e{T)V in the MCE and the 
CE. Connecting lines between the MC-dots are drawn to guide the eyes. 
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FIG. 7: The same as in Fig. IHlbut for -q meson. 
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FIG. 8: The same as in Fig. IHl but for p meson. 
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FIG. 9: The same as in Fig. IHl but for /2 meson. 
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FIG. 10: The same as in Fig. IHl but for J/il] meson. 
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FIG. 11; The energy dependence of the ratio of the mean multiplicities of (top), r] (middle) and p 
(bottom) mesons in the MCE and the CE to those in the GCE obtained within the hadron-resonance 
gas model. The vertical dashed lines indicates the threshold energies of several two particle channels 
which are responsible for the observed structures. The calculations are performed for T = 160 MeV 
and E = {E)gce in the GCE and for E = e{T)V in the MCE. Connecting lines between the MC-dots 
are drawn to guide the eyes. 
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FIG. 12; The dependence of the mean multiplicity of tt® mesons on the system energy obtained within 
the hadron-resonance gas model by use of the micro-canonical (solid line) and grand canonical (dotted 
line) ensembles. The calculations are performed for T = 160 MeV and E = {E)cce in GCE and for 
E = e{T)V in MCE. Connecting lines between the MC-dots are drawn to guide the eyes. 
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FIG. 13: The energy dependence of the ratios (^/)/(vr°), {p)/{ t^^) and {J/if))/{' k^) obtained within the 
hadron-resonance gas model by use of the micro-canonical (points and solid line) and grand canonical 
(dotted line) ensembles. The calculations are performed for T = 160 MeV and E = {E)cce in GCE 
and for E = £{T)V in MCE. Connecting lines between the MC-dots are drawn to guide the eyes. 
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